Phase Behavior of Columnar DNA Assemblies 
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The pair interaction between two stiff parallel linear DNA molecules depends not only on the dis- 
tance between their axes but on their azimuthal orientation. The positional and orientational order 
in columnar B-DNA assemblies in solution is investigated, based on the DNA-DNA electrostatic 
pair potential that takes into account DNA helical symmetry and the amount and distribution of 
adsorbed counterions. A phase diagram obtained by lattice sum calculations predicts a variety of 
positionally and azimuthally ordered phases and bundling transitions strongly depending on the 
counterion adsorption patterns. 

PACS numbers: 82.35.Rs, 64.70.-p, 87.14.Gg, 82.70.Dd 



DNA is a polyelectrolyte molecule. In aqueous elec- 
trolyte solutions, cations along its helices dissociate from 
it and dissolve into the mixture, leaving behind nega- 
tive charges that reside on the phosphates of the DNA 
backbone. According to the Manning condensation the- 
ory, a fraction of the cations condenses into the Bjerrum 
layer near the molecular surface If some of the ions 
specifically adsorb onto DNA, its surface could be almost 
fully neutralized Q or even overcharged || . Far from its 
axis, DNA can be apprehended as a charged cylinder. If 
the charge were continuously smeared, there would be 
only an electrostatic repulsion between two molecules, 
exponentially screened by the electrolyte. However, the 
net distribution of charge on the molecules is not homo- 
geneous and this can dramatically alter the interaction 
potential at intermediate distances. Indeed, in order to 
condense DNA in an aggregate, one has either to apply 
osmotic stress E) or use condensing agents, such as salts 
with Mn 2+ , Cd 2+ , spermidin, protamine or cobalt hex- 
ammine [^| cations. These cations are known to specifi- 
cally adsorb on DNA, predominantly in the DNA grooves 
||. Other counterions, such as, e.g., Ca 2+ or Mg 2+ , that 
have strong affinity to phosphates and adsorb preferen- 
tially on the strands do not induce DNA aggregation. 
Obviously, one effect of these specifically adsorbing coun- 
terions is the reduction of the net charge on the DNA. 
However, were this to be the only effect, it would have 
been hard to explain the observed sensitivity to the sort 
of counterions of DNA condensation (5| and of the meso- 
morphism of resulting aggregates @ . 

Recently a new explanation of the features of DNA ag- 
gregation was suggested [|| resting on a Debye- Bjerrum 
theory of electrostatic interaction between helical macro- 
molecules The theory offered first a formalism for a 
description of interaction between cylindrical molecules 
(with parallel axes) for arbitrary surface charge distribu- 
tions on the molecules Then it explored its conse- 
quences for helical charge distributions, including those 
typical for double stranded B- and A-forms of DNA 



H [T(J. Various patterns of adsorbed counterions, includ- 
ing those spiraling through DNA major and minor groves, 
were considered. Thus, the effect of helically structured 
separation between negative and positive charges on each 
molecule was rationalized, explaining, in particular, a 
stronger DNA-DNA attraction in the presence of coun- 
terions preferentially adsorbing into the major groove. A 
number of applications of the theory [ |ic| proved to be in 
line with experimental observations, and the main prop- 
erties of the calculated interaction potential were verified 
by computer simulations 111]]. 

Let us draw a plane perpendicular to the parallel 
axes of the molecules. For each molecule draw a vec- 
tor joining the axes where the 5' — 3' strand |L2| hits 
the plane and call it 'spin'. The angle between the 
two spins, <f), may be called the angle of mutual az- 
imuthal orientation of the two molecules. A remark- 
able effect of DNA double strandedness is a peculiar 
dependence of the interaction potential on (f>. To a 
good approximation the (^-dependent part of the po- 
tential reads g u(R,4>) = — A(R) cos 4> + B(R) cos 2 0, 
where A(R) and B(R) are positively-definite functions 
of the interaxial separation R, with A(R) dominating 
at large R. For this potential, the optimum angle is 
4> = ±Q (2B(R) ~ A(R)) cos- 1 [A(R)/2B{R)}, with the 
Heaviside step function Q(x). In other words, there are 
two symmetrical nonzero values of the angle at distances 
smaller than a critical distance at which A(R) = 2B(R), 
and zero value of the angle at large distances. The ratio 
A(R) I B(R) diminishes with decreasing R and the abso- 
lute value of the angle grows. The values of the functions 
A(R) and B(R) depend on the parameters of the DNA 
helical structure and distribution of adsorbed ions, but 
typically the absolute value of the optimum angle varies 
between and tt/2. 

Thus the problem of statistical properties of columnar 
aggregates of long DNA molecules can be mapped on a 
2d-problem of XY-spins interacting via such an unusual 
potential. Since A(R) and B(R) exponentially decay 
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with R, the dominant role is played by nearest neighbor 
interactions. While the = 0-case is compatible with 
simple a hexagonal lattice, the case <\> ^ results into 
frustrations of positional and orientational order [fl3|| . 
Due to the coupling between the positional and orien- 
tational variables in the interaction ( l R — (f> coupling'), 
one may expect most peculiar positional and spin struc- 
tures in the aggregate, a feature known as mesomorphism 
of DNA assemblies @. 

In this work, we analyze the statistical properties 
of such assemblies in aqueous solutions. We calculate 
phase diagrams that depend on the DNA- and salt- 
concentrations, and on the counterion adsorption pat- 
tern. To investigate the stability of various phases, we 
carry out lattice-sum calculations for interacting DNA 
molecules and supplement them with the cntropic and 
cohesive contributions from the ions of the solution. The 
so-obtained variational Helmholtz free energy is finally 
minimized among the candidate phases and the equilib- 
rium states are obtained. Treating DNA molecules as 
rigid is justified as long as their contour length does not 
exceed the persistence length L p — 500 A. The axes of 
the molecules remain parallel (to the z-axis) as long as 
the nearest-neighbor distance in the aggregate remains 
below 40 A. To model the interaction, we envision the 
molecules as long cylinders, carrying helical, continuous 
line charges on their surface. Each DNA-doublc helix 
carries the negative charge of phosphates with surface 
charge density a = 16.8/iC/cm 2 plus a compensating 
positive charge coming from the adsorbed counterions. 
Let < 6 < 1 be the degree of charge compensation, /i, 
/2, and fs be the fractions of condensed counterions in 
the major and the minor grooves, and on the two strands, 
respectively (fx + fa + /a = 1). The mobile counterions 
in solution screen the Coulomb interactions between the 
helices, causing at large separations an exponential decay 
of the latter with the Debye screening length Sol- 
vent screening is accounted for by its dielectric constant 
e. For DNA structural parameters we take the B-DNA 
values: pitch H ss 34 A (g = 2it/H) and the hard-core 
radius a = 9 A. For the pair interaction potential, we 
take the form B (R > 2a): 
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where Az is a vertical displacement, equivalent to a 
'spin angle' 4> = <?Az [[14). Here, uo = 8-Ka 2 /eK 2 
(« 2.9 fcsT/A at physiological ionic strength), (f> s sa 0.4 7r 
is the azimuthal half-width of the minor groove, and 
I 2 - H n , m (x,y) is given by 
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FIG. 1: A schematic view of generating candi- 
date ordered spin phases of the system. (a): for 
the .HEX-lattice; (b): for the REC- and SQ- 
lattices; (c) and (d): fir the RHO- and OBWattices. 



with the modified Bessel functions K n (x) and Ij(y). The 
primes denote derivatives. The precise form of the (jj- 
dependence is affected by the distributions /i, i — 1, 2, 3 
of the condensed counterions || . Keeping only the n = 0- 
tcrm in the sum of Eq. (Q) returns a pair potential de- 
pending on R only, as it corresponds to the approxima- 
tion of continuously charged cylinders. Truncating the 
infinite sum of Eq. (Q) at |n| = 5 is sufficient for full 
convergence of the sum, for all cases studied here. 

For all case studies of this work, the pair potential is 
greater than &bT, thus we focus on the ground state- 
analysis of the basic structures of the assembly. To this 
end, we considered the five two-dimensional Bravais lat- 
tices, i.e., the hexagonal (HEX), square (SQ), rectangu- 
lar (REC), rhombic (RHO) and oblique (OBL) lattices. 
In order to explore the ordered spin structures, we con- 
structed a certain spin pattern on the elementary pla- 
quette of every lattice and repeated it along the lattice 
directions. If all site-site interactions in the Hamiltonian 
are contained within the elementary geometrical cell (pla- 
quette) of the lattice, then the exact ground state can be 
obtained as follows. The energy of this plaquette must 
be minimized with respect to the spin angles, and then 
the optimized spin pattern on the plaquette must be re- 
peated throughout the lattice. We have interactions of 
higher-order-neighbors in our model but the exponential 
decay of the i?-dependent prefactors guarantees that the 
nearest-neighbor interaction dominates. We have kept up 
to 10 neighboring shells in the calculations of the lattice 
sums, that turn out to be sufficient for convergence. 

In Fig. ^ we show schematically the algorithms em- 
ployed for the generation of the ordered spin structures. 
Choosing the orientation of one of the spins as reference 
(4> = 0), we are left with two free orientations per pla- 
quette for the ffiJX-lattice and three for the REC- and 
(SQ-lattices. The lattice is filled by successive mirror- 
reflections of the cells across their edges, as shown in 
Figs. |l](a) and (b). As far as the RHO- and OBL- lattices 
are concerned, the procedure involving three free spin an- 
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FIG. 2: The four stable magnetic phases. The ar- 
rows indicate the relative orientations of the helical 
DNA molecules. The acronyms stand for ferromagnetic 
(FM), antiferromagnetic Ising (AFI), antiferromagnetic 
Potts (AFP), and antiferromagnetic Heisenberg (AFH). 



gles per plaquette does not generate identical plaquettes 
upon reflection: in these lattice types there is a short 
and a long diagonal which exchange their roles upon re- 
flection. We employ two complementary algorithms for 
generating ordered magnetic structures on these lattices: 
first, we place spins of orientations fa and fa along the 
cell edges and fa — fa along the long diagonal and use the 
successive reflection algorithm. This guarantees that all 
pairs of spins across all diagonals will have relative angles 
fa — fa, see Fig. |l](c). Alternatively, we place along the 
long diagonal a spin with an angle fa + fa and subse- 
quently we increase the spin angle along the horizontal 
direction by an amount of fa, and along the oblique di- 
rection by an amount of fa for every step. We generate 
thus structures in which all spins along the short diago- 
nals have an angle fa — fa and all spins along the long 
diagonals an angle fa + fa, as shown in Fig. |l|(d). 

The interaction between any two molecules is 
L p u(R,faj, where L p = 500 A is the persistence length 
and u(R,fa) is given by Eq. (Q). It was found that the 
energy needed to destroy the translational or orienta- 
tional order must be more than several k B T at room 
temperature, hence the lattice-sum calculations provide 
the representative thermodynamic states. The 2d DNA- 
concentration p was varied within < pa 2 < l/(2\/3), 
the upper limit corresponding to the close-packed config- 
uration in a HEX-lsXtice. For every density, minimiza- 
tions of the lattice energy with respect to the plaquette 
sets {fa}, the size ratios b/c (for the i?£'C-lattice) and/or 
the geometrical angle u> (RHO- and 05.Mattices) , Fig. 
|l|, were carried out. This way the optimized lattice-sum 
energy, Ux(&,p), was obtained, where X stands for the 
lattice type and €> = (fa, fa, . . . , faj) denotes the config- 
uration of the N spins in the system. 

To access the full thermodynamics of the DNA 
solution-salt mixture, we have to add the contributions 
to the free energy arising from the counter- and co-ions, 



(numbers N± and concentrations c±, respectively.) The 
effect of these degrees of freedom is to add an exten- 
sive term to the free energy of the system Jlj|, F c — 
F° + F° + F coh , where F£ = N± k B T [ln(c±A|) - l] 
are entropic contributions from the kinetic part of the 
Hamiltonian with the thermal de Broglie wavelengths A± 
of the counter- and co-ions, and 
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is a cohesive term. In Eq. (||), e is the electron charge, 
Z\e\ — 2iraL p a(l — 9) is the uncompensated DNA-charge, 
c + = Zp/L p + n s and c_ — n s , with the salt concen- 
tration n s . Finally, V is the volume of the system and 
k = ^4tt(Z pj L p + 2n s )e 2 / (ek B T) for monovalent salt 
ions. The Helmholtz free energy is Ux + F c . 

When counterions are condensed on strands, i.e., f\ = 
f-2 = and / 3 = 1, the DNA-DNA interaction is purely 
repulsive. The system is found to crystallize into the 
HEX lattice at all DNA-densities but a large variety of 
orientational (magnetic) structures occur, as a result of 
the frustration of the system. The structures are shown 
in Fig. U and the phase diagram of the DNA-salt mix- 
ture in Fig. ||(a). The phase denoted FM is a simple 
ferromagnetic phase, in which all DNA-molecules have 
the same azimuthal orientation. The phase denoted AFI 
displays antiferromagnetic-Ising type ordering, with half 
DNA-molecules having a given orientation angle if they 
lie on one of the sublattices and a different orientation 
on the other. The AFP phase has a three-state antifer- 
romagnetic Potts |ll| type of ordering, with 1/3 of the 
spins pointing in a reference direction = 0, 1/3 in the 
angle fa and 1/3 in the angle 2fa. Note that the angle fa 
grows with DNA concentration. Finally, the AFH-pha.se 
has the orientational ordering of the two-dimensional an- 
tiferromagnetic Heisenberg model, with spins residing in 
the three sublattices of the hexagonal lattice having mu- 
tual orientational angles of 120° to one another. The 
transition between the FM and AFP phases is second- 
order but the AFP — > AFI and AFI — > AFH transitions 
are first order. Referring to Fig. ||(a), we see that the 
FM phase is stable at low DNA-concentrations. Indeed, 
for such average intermolecular separations the optimal 
orientation angle between the molecules is zero. The non- 
trivial phases arise at higher densities of the aggregates. 
Similar mcsophases were found recently within the frame- 



work of a phenomenological Landau theory |17 



When counterions condense in grooves, an attraction 
between the DNA-molecules arises, since the possibility 
of having positively charged parts of one molecule ap- 
proaching close to negatively charged parts of the other 
through an appropriate mutual orientation opens up. 
This leads to broad phase coexistence lines between dense 
DNA-aggregates and DNA-free solutions. This is demon- 
strated in Fig. |^(b) for the case f\ — 0.3, f% = 0.7 and 
/.3 = for 9 — 0.9. Lowering 9, i.e., increasing the 
Coulomb repulsion, opens up a hexagonal phase at small 
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FIG. 3: Phase diagrams of DNA-salt mixtures obtained by the procedure described in the text: (a) 8 — 0.9, fs = 
1. The geometrical lattice is here HEX. (b) 6 = 0.9, /i = 0.3, / 2 = 0.7; (c) 6 = 0.7, /i = 0.3, f 2 = 0.7. 
Dashed lines denote second-order magnetic transitions and solid lines first-order ones. The geometrical transitions 
between different lattice types in (b) and (c) are 2nd order; the straight lines are tielines between coexisting phases. 



DNA- and electrolyte-concentrations, see Fig. ||(c). In 
the one-phase region, a rhombic phase shows up for mod- 
erate to high densities and, due to packing constraints, a 
HEX crystal appears at very high DNA-concentrations. 
A strong qualitative difference in the macroscopic be- 
havior of columnar DNA assemblies arises, depending 
on whether the countcrions condense on strands or in 
grooves. In the former case, all transitions are 'magnetic' 
in nature. In the latter, DNA-bundling takes place and 
rhombic lattices are stabilized. 

The predictions of the theory ask for experimental ver- 
ification. Such a task is not easy, since the reliable data 
up to date refer only to highly concentrated phases [ pj| , 
where the number of the basic assumptions inherent to 
the form of the pair potential may be questioned (the 
Dcbyc-Bjcrrum approximation, independence of solvent 



dielectric constant on the aggregate density, effects of 
nonlocal polarizabilty, etc.). Such a verification might 
become possible with the increase of experimental res- 
olution in X-ray diffraction, which could open the way 
for the study of less dense aggregates. The predicted 
specific effect of cation adsorption on the phase dia- 
gram is particularly challenging. Since the adsorption 
isotherms and the distributions of the adsorbed ions be- 
tween the minor and major grooves are poorly known, 
one should concentrate here on the qualitative effects, 
i.e., the (dis) appearance of mesophases triggered by the 
presence of different DNA condensing counterions. 
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